The decay law can have an irregular character 
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Within a well-known decay model describing a particle confined initially within a spherical 5 po- 
tential shell, we consider the situation when the undecayed state has an unusual energy distribution 
decaying slowly as fc ^ oo; the simplest example corresponds to a wave function constant within 
the shell. We show that the non-decay probability as a function of time behaves then in a highly 
irregular, most likely fractal way. 
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The decay of an unstable quantum systems is one of 
the effects frequently discussed and various aspects of 
such processes were considered. To name just a few, re- 
call the long-time deviation from the exponential decay 
lawQ, the short-time behavior related to the Zeno and 
anti-Zeno effects HjSE^ revival effects such the clas- 
sical one in the kaon-antikaon system, etc. In all the 
existing literature Q, however, the decay law is treated 
as a smooth function, either explicitly or implicitly, e.g. 
by dealing with its derivatives. The aim of this letter is 
to show there are situation when it is not the case. 

A hint why it could be so comes from the behavior 
of Schrodinger wave functions during the time evolution. 
While in most cases the evolution causes smoothing 8], 
it may not be true for for a particle confined in a po- 
tential well and the initial state does not belong to the 
domain of the Hamiltonian. A simple and striking ex- 
ample was found by M. Berry 9] for a rect ang ular hard- 
wall box, and independently by B. Thaller for a one- 
dimensional infinite potential well. It appears that if the 
initial wave function is constant, it evolves into a steplike- 
shaped iplx^t) for times which are rational multiples of 
the period, t = qT with q = N/M, and the number of 
steps increases with growing M, while for an irrational q 
the function ')p{x, t) is fractal w.r.t. the variable x. 

One can naturally ask what will happen if the hard wall 
is replaced by a semitransparent barrier through which 
the particle can tunnel into the outside space. In a broad 
sense this is one of the most classical decay model which 
can be traced back to 0. We will deal with its par- 
ticular case when the barrier is given by a spherical 5 
potential which is sometimes called Winter model b eing 
introduced for the first time, to our knowledge, in 
see also |0 . The described behavior of the wave function 
in the absence of tunneling suggests that in the decaying 
system the irregular time dependence could also be visi- 
ble, both in the wave function and in various quantities 
derived from it at least in the weak coupling case. 
The aim of the present letter is to demonstrate that this 
conjecture is indeed valid. 

To be concrete, we will study a spinless nonrelativistic 
quantum particle described by the Hamiltonian 



with a fixed i? > 0; we use rational units, h — 2m — 
1. For simplicity we restrict our attention to the s-wave 
part of the problem, writing thus the wave functions as 
ip{r,t) — -^=r~^(j){r, t) with the associated Hamiltonian 



ha = --r^ + oi5(r ~ R) 



(2) 



in the lowest partial wave. We are interested in the time 
evolution determined by the Hamiltonian ij}{f,t) = 
e~'^°*V'(^i 0) for a fixed initial condition il){r, 0) with the 
support inside the ball of radius R; the advantage of the 
used model is that the propagator can be computed ex- 
plicitly. Of a particular interest is the decay law. 



P{t)^ / |0(r, t)|2dr, 



(3) 



i.e. the probability that the system localized initially 
within the shell will be still found there at the measure- 
ment performed at an instant t. We are going to derive 
an exact formula for the decay law which we will then 
allow us to evaluate the function Q numerically for a 
given initial state. 

It is straightforward to check |lj] that the Hamilto- 
nian iQJ has no bound states. On the other hand, it 
has infinitely many resonances with the widths increas- 
ing logarithmically w.r.t. the resonance index. A natural 
and well known idea is to employ them as a tool 

to expand the quantities of interest. 

First of all, we have to find Green's function g{k, r, r'), 
i.e. the integral kernel of (h^ — A;^)^^ which determines 
the time evolution in the standard way [l7j . 
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recall that u{ha) = [0,oo) for a > 0; we are going to 
perform a resonance expansion of the integral The 
Green function for a system with singular potential is 
obtained from Krein's formula. 
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where is the free Green function, in particular, 

<i)^(r) = i sin(fcr) e''^''^ holds for r < R, and A(fc) is de- 
termined by boundary conditions at the singular point 
i?; by a direct calculation |0 one finds 



m = - 



2k 



[1 - e^'^R) 



(5) 



This allows us to write the sought reduced wave function 

oo 

at time t as 0(r, t) — J u{t, r, r')(f){r' , 0) d?-' with 



u(t, r, r ) — — lim e lmg(k + ie, r, r ) 2k dk 
TT eio J 


Since has no eigenvalues we may pass in the last for- 
mula to the limit e — s- obtaining 



u{t, r, r') = / p{k, r, r') exp{-ik'^t) 2k dk, (6) 



where p{k, r, r') = ^Im g{k, r, r'). This can be using 
equation ((SJ written explicitly as 



p{k, r, r') 



2k sin(fcr) sin(fcr') 



7r(2fc2 + 2a2 gin^ kR + 2ka sin 2kR) 



(7) 



The resonances of the problem are identified with the 
poles of g{-, r, r') continued analytically to the lower 
momentum halfplane. They exist in pairs, those in the 
fourth quadrant, denoted as fc„ in the increasing order of 
their real parts, and — fc„. The set of singularities of the 
kernel then include also the mirrored points, being 



S — {kji, kn, fc„, kji : ri G N}. 



(8) 



In the vicinity of the singular point fc„ the function 
p{-, r, r') can be written as 

/, /N * Vnir)vn{r') 

r, r ) = — — p hx(fc, r, r), (9) 

where w„(r) solves the differential equation haV„{r) = 
k'^Vn(r) and the function x is locally analytic. 

The factor ^ is chosen to get the conventional normal- 
ization of the resonant state Vn{r) [l5| . The last named 
paper also demonstrates that for a finite-range potential 
barrier and r, r' < R the function p{-, r, r') decreases in 
every direction of the fc-plane; in the present case it is 
not difhcult to verify this claim directly. Then one can 
express p(fc, r, r') as the sum over the pole singularities 



p(fc, r, r') 



1 



fees 



k-k 



:.Res^p(fc, r, r') (10) 



and derive from the residue theorem the following useful 
formula 
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ReSj,p(fc, r, r') = . 



(11) 



The last two equations can be rewritten in view of eq. 
and the symmetry of the set S in the form 

p{k, r, ^0 = E ^ fc2 ^ fc2 ^ Vn{r)vn{r') , (12) 



ne2 



Y,^^n{r)vn{r')^0, (13) 



where we denote fc_„ := —kn and u_„ is the associated 
solution of the equation HQU_„(r) = fc^„u_„(r). 

Next we substitute from eq. p2(l into JHl and using 
the identity H13|l we arrive at the formula 



u{t, r, r') = 



E 



exp{-ik'^t) 2fc2 



Vn{r)vn{r')dk 
M{kn, t)vjr)vn{r') (14) 



StT / ^ k-ri 

riSZ " 



with M{kn, t) = ^ e"" erfc(M„) and u„ := 



"/^k^Vi. 



Indeed, using 2k'^ = 2{k'^ — fc^) -I- 2kn we write u{t, r, r') 
as a sum of two terms the first of which vanishes in view 
of H13|) . The second one decomposes again into a sum 
of two integrals containing fc„ ± A; in their denominators, 
which gives the right-hand side of H14|) with 



I k kfi 



■dk. 



(15) 



in other words, the above expression. 

Now a straightforward calculation using Q and H14|l 
allows us to express the decay law in the form 

P{t) = CnClInlM{k^, t)M{ki, t) , (16) 



with the coefficients 

R R 
C-a ■■= J Hr, 0)v„(r)dr, := J i;„(r)w, (r) dr . (17) 



This expression holds generally, see [T^ : in our particu- 
lar case of the Hamiltonian (0) we can specify w„(r) — 
V2Qn sin(fc„r) with the coefficient Qn equal to 



'2iki 



2kn + a^i? sin 2fc„i? -|- a sin 2knR + 2fc„ai? cos 2fc„ J? 



1/2 



Now we are ready to compute the decay law for a 
given initial state. Without loss of generality we may 
put R = 1, we choose the value a = 500 for numerical 
evaluation and replace the infinite series by a cut-off one 
with \n\ < 1000. As the first example we consider initial 
wave function constant within the well, i.e. 



{r, 0) = R-^/^V3r , r<R. 



(18) 
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FIG. 1; Decay law for the initial state (j>{r, 0) = R'^^^VSr. 
In the inset we plot the logarithmic derivative averaged over 
intervals of the length approximately r/200. 
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FIG. 2: Decay law for initial state c/>(r, 0) = R and its 
logarithmic derivative, locally smeared, in the inset. 



The corresponding decay law is plotted in Figure Q It 
is irregular having numerous steps the most pro- 
nounced at the period T = 2R^ /tt and its simple ratio- 
nal multiples. To make them more visible we plot in the 
inset the logarithmic derivative of the function P{t); it is 
locally smeared, otherwise the picture would be a fuzzy 
band. The irregular structure is expected to be fractal; 
it persists at higher time but its amplitude decreases rel- 
atively w.r.t. the smooth background. 

In the next example we choose the initial state hav- 
ing constant reduced wave function, (j){r, 0) = i?^^/^ for 
r < R, so ^{r^ 0) has a (square integrable) singularity at 
the origin; the advantage is that the reduced problem of- 
fers a straightforward comparison to the one-dimensional 
example treated in [lol | including the shapes of the wave 
functions. The decay law for this case is plotted in Fig- 
ure [3 it again exhibits derivative jumps around simple 
rational multiples of the period. The corresponding func- 



FIG. 3: The probability density inside the sphere with R = 1 
muhiplied by at the instants t = T/S, r/16, and T/27. 



tion |(/)(r, t)p for three such values is plotted in Figure 13 
To compare with notice that the one-dimensional 
well has length i? = 1 so the revival period of this state 
is r/8. In the absence of decay the function for t = T/8 
is just constant, the other two are simple step functions. 
We see that the tunneling through 5-barrier modifies the 
shape of the function mostly in the vicinity of the origin, 
the barrier, and the jump points. 

To support the conjecture about the fractal character 
of P{t) let us look how its derivative behaves in the limit 
a ~* oo when (j)(r, t) expands asymptotically as 



(r, t) 



|C„exp(-ifc^t)u„(r). 



(19) 



It is easy to see that for a fixed n and a —f oo the reso- 
nance position expands around fc„^o nir/R as 



kn,0 77 

an 



(ai?)2 'a^i? 



(20) 



which shows that in the leading order we have Vn{r) ~ 
sin(fc„r), and furthermore, that the substantial con- 
tribution to the sum in H19|l comes from terms with 
n < [a^~^^] for some < e < 1/3. 

The derivative P{t) can be computed from the proba- 
bility current conservation _7(r) = integrat- 
ing it over the interval (0, R) and using j(0) = we get 



P{t) = -2Im ((/)'(i?, t)(j){R, t)) . (21) 
We plug the above expansion into 119|l obtaining 



cp{R, t) 



'-^(-l)"C„exp 



exp 



9P 



a oP- 



(22) 
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and a similar expansion for (t>'{R, t) with the last bracket 
replaced by fc„^o- We observe that for j > — 1 we have 

where we have denoted Ij :— J e~^' Ax — {'^-y')- 



Let us now we assume that the coefficients in (|19|l sat- 
isfy Cn ^ fc^o as n ^ oo. Suppose first that the decay 
is fast enough, p > 1; notice that this is certainly true in 
the finite-energy case with p > 3/2. The term —knfi/a 
in H22|) obviously does not contribute to the imaginary 
part, hence we find that \P{t)\ < const a4/3-4/3p _^ q 
holds as a — > oo uniformly in the time variable. 

The situation is different if the decay is slow, p < 1. 

As an illustration take C„ = (— 1)"^^"^^"' "^hich corre- 
sponds to the first one of the above numerical examples. 
Since the real part of the resonances changes with a, 
cf. H2UI) . it is natural to study the limit oiP{ta) as a — *■ oo 
at the moving time value ta := t{l + 2/aR). Up to higher 
order terms the appropriate value (p^R, ta) is obtained by 
removing the bracket (1 — 2/aR) at the right-hand side 
of H22() and (j)'{R, ta) is obtained similarly. 

Consider first irrational multiples of T. We use the 
observation made in j^] that the modulus of X]n=i e*'^"^* 
is for an irrational t bound by C L^~^ where C, s depend 
on t only. In combination with a Cauchy-like estimate, 
"«^n < Y,'^=i I J2]=i "il l^n - bn+i\ which yiclds 

5^exp(-zfc^ ot)exp ^t ki^, < const a^/^W+i"-) 



and consequently, P{ta) — > as a oo similarly as in 
the case of fast decaying coefficients. 

Let us assume next rational times, t = ^T. li pq is 
odd then SL{t) J2n=i e""^* repeatedly retraces by 
[loj l the same pattern, hence P{ta) — > - cf. Fig. ^ at 
the half period. On the other hand, for pq even |S'l(OI 
grows linearly with L, and consequently, lim P{ta) > 0. 

a — *oo 

As a example let us compute this limit for the period T, 
i.e. p = q ^ 1. Using H21|) we find 

lim P{Ta) - lim Im ( V exp ( ) 

= -f (f) 217^^^^° = 

notice that without a local smearing P{T) in the inset of 
Fig. ^ would take approximately this value. 

In conclusion, we have reexamined time decay in Win- 
ter model and found indications that the decay law is a 
highly irregular function if the energy distribution decays 
slowly as fc ^ oo. 
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